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Abstract 

We define generalized bivariate polynomials, from which specifying 
initial conditions the bivariate Fibonacci and Lucas polynomials are 
obtained. Using essentially a matrix approach we derive identities and 
inequalities that in most cases generalize known results. 


1 Antefacts 

The generalized bivariate Fibonacci polynomial may be defined as 
H n (x, y ) = xH n _ i(x, y) + yH n _ 2 {x, y), H 0 (x, y) = a 0 , Hi(x, y) = cp. 

We assume y ^ 0, as well as x 2 + 4y ^ 0. 

If we set ao = 0, ai = 1 we obtain the bivariate Fibonacci polynomials 
F n (x, y)\ with ao = 2, a\ = x we get the bivariate Lucas polynomials. 
L n (x, y). The univariate polynomials are obtained setting y = 1. The 
characteristic polynomial is 


fit) = t 2 -tx-y, 

and the roots are (see for example [3]) 

x + y x 2 + 4y oo, \ x ~ \/x 2 + 4y 
a = a{x,y) = ---, P =/3{x, y) = ---• 

Note that a + (3 = x, a(3 = —y. Furthermore a > /3; if y > 0 then a > 0 
and (3 < 0 and if x > 0 then a > \(3\, while if x < 0 then a < \(3\\ if y < 0 
then a and /3 have the same sign: if x > 0 they are both positive, while if 
x < 0 they are both negative. 
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The ordinary generating function g(t ) is 


ait) 


a 0 + («i - a 0 x)t 
1 — xt — yt 2 


which can be obtained through the usual trick: write 


g(t ) — Hq + H\t + H 2 F + • • • 


then obtain xtg(t) and yt 2 g{t ), subtract them from g(t ), use the recurrence to 
get rid of all summands except the first two and impose the initial conditions. 
Using the Rational Expansion Theorem (see [5]) we get the Binet’s form 


H n (x, y) 


(ai - /3a 0 ) 


a 


n 


a — (3 


((il — Clflo) 


(3 n 

a — (3 


(1) 


Hence 


F n (x, y) = 


a 


(3 r 


a — /3 a — P 
L n (x, y) = a” + (3 n . 


In the sequel we will write simply H n , F n and L n instead of H n (x, y), F n (x, y) 
and L n (x , y). 

The sequence can be extended to negative subscripts by defining 


H 


—n 


- H-(n- 1) H-^-(n-2)- 

y y 


In this case the characteristic polynomial is the reflected polynomial (see 
[5]) of fit): hence the roots are ^ and jj. Define the matrix 


A = 


x 1 

y o 


Note that |A| = —y. 

By easy induction 

A n 


F n +l F n 

yp'n y Fn i 


From the above equation taking determinants 


F n ~iF n +i — F 2 


(-I)V -1 - 


( 2 ) 


(3) 

(4) 
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Also from A 2n = A n A n it follows 

Fn +1 + yFn = ^2n+l- (5) 

The eigenvalues of A are a and (3: from this 

tr(A n ) = F n+1 + yF n . _i = a n + f3 n = L n . (6) 

Also it is immediate to get 

FnF n F 2r 1- 

A wealth of inequalities can be obtained from Equation 3 using the classical 
matrix inequality on the extremes of quadratic forms 

f3 n (zl + z\) < z\F n+l + ziz 2 {l + y)F n + < a n (z{ + z$), (7) 

for all real z i, z 2 . 

Now define the matrix 


a\x + a 0 y a i 
a\V a 0 y 
= a 0 yI + aiA. 

Again by easy induction we get 

Hn+2 H n +1 
u F n ■ i yH n 

Taking ao = 2, ai = x we get the matrix 

B _ x 1 + 2y x 
xy 2 y 
= 2yl + x’A, 




and 


Because 


_ -^n+2 i'n+l 

yFn+i yF n 

CA n = (o 0 yl + aiA)A n = a 0 yA n + aiA n+1 , 


( 8 ) 
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we obtain the identity 


Hji+i — oq yF n Qj\F n +h 


which specializes to 

L‘n +1 = ^yF'n + xF n + 1 . 

Analogous results we get for //_„ using the inverse of A 


A' 1 = 


0 


l 

y 

1 -2 

y 


In particular we get 


(-l) n 

H—n = [flO-Pn+l ®l-^n] • 

y 


Because |C| = aociixy + a^y 2 — a\y taking determinants in Equation 8 we 
get 

H n H n+2 - H 2 +1 = (-l) n y n ^ 1 (a 0 aixy + a^y 2 - a\y ), 

which can be termed the Simpson formula. 

Using identity 1.60 in [4], that is 


LfJ 

£(- 1 ) 

k =0 


, n — h\ ™n+l _ n +1 

v #)”-“ = -- 


writing a for x and /3 for y we have 

LfJ 


£(-i)* 

k=0 


n — k' 
, k i 


(cx(S) '(a + /3) 


n—2k 


x-y 


n +1 _ / 3 n+l 


« n+i - 0 
a — (3 

= F n + 1, 


( 10 ) 


which means 


LfJ 

£ 

k=0 


n ~ k k )y k x n - 2k = F n+1 . 
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2 Inversion formulas 

Now we will establish the following formula for 7 complex and i 2 = — 1 


Fn(x, -y) = (ry) (n l) F n [i^x, 7 2 y). 


We have 


where 


F n (x, -y) = 


a 

a\ - Pi 


x + 1/x 2 — 4 y 


Analogously for P\. Also 


where 


tp r 2 \ a 2 ~ P 2 
F n {yyx, 72/) =-7-, 

«2 - P 2 


a 2 = 


i'fx + y —7 2 x 2 + 47 2 y 
2 

i7x + i'j\/x 2 — 4 y 


. x + \Jx 2 — Ay 

= 17-j- 

= i7Q]. 


Analogously for P 2 . Then 


F n (Fj/x, 7 2 y) = 


(i7) n < - (i7) n ,Sf 


* /t \ / 7 / £7 / • • O 

17Q!l — 17Pl 

= (i 7 ) n ^ 1 F n (x, -y), 

whence the desired formula. Some consequences. If we set 7 = i we get 


F n (x, -y) = (-l) n 1 F n (-x, -y ). 


Also from 


F n {x, y) = (hy) (n 1 ] F n (i^x, -7 2 y) 
we get with 7 = 7 = = 


F r) (x,y) = (Vy) n - i T ri —,1 
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In a similar fashion we obtain 


and 


It follows 


Ln(x, -y) = ( 17 ) n L n { vyx, 7 2 y), 
L n (x 1 y ) = (Vy) n L n lj , 


L n (—x, y) = (—l) n L n (>, y). 


(14) 


F 2n+ i(-x, y) = F 2n +i{x, y), F 2n (-x, y) = -F 2n (x, y), 


L 2n +i(-x, y) = -L 2n +i{x, y), L 2n (-x, y) = F 2n (x , y). 

From the Binet’s forms we get 


F n (ax, by) = a n l F n 



from which we get 


L n (ax , by) = a n L 




F n (x, y) = x n l F n ^1, , 

L n (x, y) = x n L n ^1, . 


3 A complex representation 


Consider the matrix 


B = 


\Jx + 2i 


x + i 1 
1 i 


where i = >/—1. Then 


This matrix is obtained by Equation 2 setting y = 1 and so is the matrix 
generating the univariate Fibonacci polynomials. Then 


F n+ i(x, 1) F n (x, 1) 
F n (x, 1) F n _ i(x, 1) 
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Using the method exposed in [7] to B 2n = A n (which consists in expressing 
the elements of the n-th power of a 2 x 2 matrix in terms of the elements of 
the matrix) we get 

n ~ 1 (‘In — 1 — r\ 

F„(z,l) = £ + (15) 

r=0\ 7 / 

Using Equation 13, after replacing in the previous equation x by we get 
n ~ l (2n — — r\ 

F n (x, y) = £ (x + 2^) n - r ~\-Vyi) r - (16) 

r=0 V ? / 


4 Summations 


Since 



k =o 


1 - a n+1 
1 — a 


(analogously for (3) and (1 — a)(l — /3) = 1 — x — y, we have 




fc =0 


ai - Qq/3 
a — (3 



k =0 


ai — aoa 

a — (3 




k =0 




a 0 (l - x) + ai - yH n - ffn-i 
1 — x — y 


which can be rewritten as 


(1 - x - y)^2 H k = a 0 (l - x) + a x - yH n - H n -i. 
k =o 


Specializing we get 


(1 - x - y) E F k = 1 - yF n - F n+1 , 
k =0 


Since 


(1 ~ x — y) E Lk = 2 — x — yL n - L„+i. 
fc =0 


E« 2fc 

k=0 


1 - a 2n+2 
1 — a 2 ’ 
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(analogously for (3) , (1 — a 2 )(l — (5 2 ) = 1—x 2 —2 y+y 2 , and a 2 —/3 2 = x(a—(3) 
we have 

n 

(1-2 y- x 2 + y 2 ) ^ F 2fc = x + y 2 F 2n - ^+ 2 , 


(1 - 2 y - x 2 + y 2 ) Y L 2k = 2 - x 2 - 2y + y 2 L 2 n - -^ 2 n+ 2 - 
k =0 

Just a little bit more involved is the following inequality. Consider 

L™ = (a n + (3 n ) m 

m / \ 

= E I 


Now upon expanding the summation and grouping the first and the last ele¬ 
ment, the second and the next-to-last element and so on we get the following 
formulas depending on the parity of m 

Lt.J 

L n =Y y k J(-y) knL n(m- 2 k) tf m ™ odd, 


L n = Y( : ) (- y) kn L n{ m- 2 k ) + ( Z) ("V)^ * 


m is even. 


If we define 


El = [yj + 1 [-1 + (-l) m+1 ] , £2 = 1 [(-l) m + 1] , 

we can group the two identities into the following one 


K = Y,[":^-y) knL nim-2k ) + e 2 r (-y)T 


ErJ^=4, 




Other consequences are 



Indeed, for instance the second one, we have 



E y k j * k y n ~ k <* k + E [ k ) xk y n ~ k P 

(y + xa) n + {y + xf3) n , 


and y + xa = a 2 , y + xf3 = /3 2 . 


5 Eigenvectors and eigenvalues 

Now V k, n 

A k BA n = A k (yI + A 2 )A n 
= yA n+k + A n+k+2 
= (yl + A 2 )A n+fe 
= BA n+fc . 

This means that Vfe, n the matrices and BA" commute: then (see [1], 
[2], [8]) the matrices possess the same eigenvectors. 

The eigenvectors (not normalized) of A are 

a (3_ 

zi = y , and Z 2 = y 

1 J ’ [ 1 

For the eigenvalues of BA” since 

BA” = yA n + A”+ 2 

we have 

A, = ya n + a n+2 , 

A 2 = y( 3 n + ( 3 n+2 . 

Consider (BA”) 2 = BA”BA”. Since A k and BA” commute, this means 
that we can write (BA”) 2 = BBA 2 ”, that is 

2 y + x 2 x L 2n +2 L 2n +i 

xy 2 y yL 2n + 1 yL 2n 

Using repeatedly the recurrence for L n we get 

Ln +2 + y^n +1 = ^2n+4 + yL 2n +2- (17) 
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Another immediate consequence is the following. Since 


Fn+iLn +2 T yF n L n +1 F n +\L n - |_i + yF n L n 
yF n Ln +2 A y F n —\L n +i yF n Ln+ 1 “i - y F n —iL n 

and A'"BA n = BA 2n it follows equating the elements in position (2, 2) 

L 2 n — F' n Ln +1 T yFn—\L n . (18) 

Because of the Schur inequality (see [1], [6]), that is: the sum of the squared 
elements is greater or equal to the sum of the squared eigenvalues, we have 
the following inequality 

Fl +l + (1 + y 2 )F 2 + y 2 F 2 _, > a 2n + (3 2n 

= L 2n . 

Note that the equality sign holds if y = 1. 

In the same way we get 

y 2 L 2 + (l + y 2 )L 2 +1 + L 2 +2 > (ya n + a n+2 ) 2 + (y/T + /T+ 2 ) 2 

= y~L 2n + 2yA2n+2 + L 2 n+4i 

and the equality sign holds if y = 1. 

Consider 

B 2 A” = y 2 A n + A n+4 + 2yA n+2 
= yBA n + BA n+2 

yL n + 2 + L n+ 4 yL n + i+ L n+3 

y 2 L n + 1 + yL n+ 3 y 2 L n + yLn +2 

Taking the determinants we get 

(-l) n y n (x 2 y + 4y 2 ) 2 = y 3 L n L n+2 + y 2 L 2 l+2 + y 2 L n L n+ 4 + yL n+2 L n+A 

~y L n + 1 — y^ j n +? J ~ 2y T n+ iL n+3 . 

Using repeatedly identity 10 we obtain 

A n L n+4 - T n+1 L n+3 = (-l) n (*V + 5*V +1 + 4y” +2 ). (19) 



A"BA" = 
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